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Superconductive  quantum  circuits  comprise  quantized  energy  levels  that  may  be  coupled  via 
microwave  electromagnetic  helds.  Described  in  this  way,  one  may  draw  a  close  analogy  to  atoms 
with  internal  (electronic)  levels  coupled  by  laser  light  fields.  In  this  Letter,  we  present  a  superconductive 
analog  to  electromagnetically  induced  transparency  that  utilizes  superconductive  quantum  circuit 
designs  of  present  day  experimental  consideration.  We  discuss  how  a  superconductive  analog  to 
electromagnetically  induced  transparency  can  be  used  to  establish  macroscopic  coherence  in  such 
systems  and,  thereby,  be  utilized  as  a  sensitive  probe  of  decoherence. 
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Superconductive  quantum  circuits  (SQCs)  comprising 
mesoscopic  Josephson  junctions  can  exhibit  quantum  co¬ 
herence  among  their  macroscopically  large  degrees  of 
freedom  [1].  They  feature  quantized  flux  and/or  charge 
states  depending  on  their  fabrication  parameters,  and  the 
resultant  quantized  energy  levels  are  analogous  to  the 
quantized  internal  levels  of  an  atom.  Spectroscopy,  Rabi 
oscillation,  and  Ramsey  interferometry  experiments  have 
demonstrated  that  SQCs  behave  as  “artificial  atoms” 
under  carefully  controlled  conditions  [2-8].  This  Letter 
extends  the  SQC-atom  analogy  to  another  quantum  opti¬ 
cal  effect  associated  with  atoms:  electromagnetically 
induced  transparency  (EIT)  [9,10].  We  propose  the  dem¬ 
onstration  of  microwave  transparency  using  a  supercon¬ 
ductive  analog  to  EIT  (denoted  S-EIT)  in  a  super¬ 
conductive  circuit  exhibiting  two  metastable  states  (e.g., 
a  qubit)  and  a  third,  shorter-lived  state  (e.g.,  the  readout 
state).  We  show  that  driving  coherent  microwave  transi¬ 
tions  between  the  qubit  states  and  the  readout  state  is  a 
demonstration  of  S-EIT.  We  further  propose  a  means  to 
use  S-EIT  to  probe  sensitively  the  qubit  decoherence  rate; 
the  philosophy  is  similar  to  proposed  EIT-based  mea¬ 
surements  of  phase  diffusion  in  atomic  Bose-Einstein 
condensates  [11]. 

The  three-level  A  system  illustrated  in  Eig.  1(a)  is  a 
standard  energy  level  configuration  utilized  in  EIT  [9,10]. 
It  comprises  two  metastable  states  |1)  and  |2),  each  of 
which  may  be  coupled  to  a  third  excited  state  |3).  In 
atoms,  the  metastable  states  are  typically  hyperflne  or 
Zeeman  levels,  while  state  |3)  is  an  excited  electronic 
state  that  may  spontaneously  decay  at  a  relatively  fast  rate 
r3.  In  an  atomic  EIT  scheme,  a  resonant  “probe”  laser 
couples  the  |l)<->  |3)  transition,  and  a  “control”  laser 
couples  the  |2)  <-►  |3)  transition.  The  transition  coupling 
strengths  are  characterized  by  their  Rabi  frequencies 
0^3  =  “dys  ■  Ej3  for  j  =  1,  2,  respectively,  where  d^3 
are  the  dipole  matrix  elements  and  £^3  are  the  slowly 


varying  envelopes  of  the  electric  fields.  Eor  particular 
Rabi  frequencies  fljj,  the  probe  and  control  fields  are 
effectively  decoupled  from  the  atoms  by  a  destructive 
quantum  interference  of  the  two  driven  transitions.  The 
result  is  probe  and  control  held  transparency  [9,10].  In 
recent  experiments,  ultraslow  light  propagation  due  to 
EIT-based  refractive-index  modifications  in  atomic 
clouds  has  been  observed  [15-17]. 


(a)  1 3)  (b)  DC  SQUID 


FIG.  1.  (a)  Energy  level  diagram  of  a  three-level  A  system. 

EIT  can  occur  in  atoms  possessing  two  long-lived  states  1 1),  |2), 
each  of  which  is  coupled  via  resonant  laser  light  fields  to  a 
radiatively  decaying  state  |3).  (b)  Circuit  schematic  of  the 
persistent-current  (PC)  qubit  and  its  readout  SQUID, 
(c)  One-dimensional  double-well  potential  and  energy-level 
diagram  for  a  three-level  SQC.  Using  PC  qubit  parameters 
[12-14],  we  calculate  W2  —  coi  =  (Itt)  36  GHz  and  coj  — 
CO2  =  (27r)  32  GHz.  The  simulated  matrix  elements 
(pi  sin(277/  +  2(f),fl\q)  for  (p,  q)  =  (1,  2),  (2,  3),  and  (1,  3)  are, 
respectively,  0.0704,  —0.125,  and  0.0158. 
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SQCs  also  exhibit  A-like  energy  level  structures 
[7,12,13,18-21].  One  example  is  the  PC  qubit,  a  super¬ 
conductive  loop  interrupted  by  two  Josephson  junctions 
of  equal  size  and  a  third  junction  scaled  smaller  in  area  by 
the  factor  a  <  \  [Pig.  1(b)]  [14,22].  Its  dynamics  are 
described  by  the  Hamiltonian 

=  [^|  +  (l  +  2a)^c>^] 

+  £'y[2-l-  a  —  2cos(ppCOS(p^  —  acos(27r/-l-  2^^)], 

(1) 

in  which  C  is  the  capacitance  of  the  larger  junctions, 
‘Pp.m  =  {‘Pi  —  ‘P'l)I%  <Pi  is  the  gauge-invariant  phase 
across  the  larger  junctions  i  =  {1,2},  Ey  is  the  Joseph¬ 
son  coupling  energy,  and  /  is  the  magnetic  flux  through 
the  loop  in  units  of  the  flux  quantum  Oq  [14].  Near  /  ~ 
1/2,  the  qubit  potential  landscape  [second  term  in  Eq.  (1)] 
assumes  a  double -well  profile.  Each  well  corresponds  to  a 
distinct  classical  state  of  the  electric  current,  i.e.,  left  or 
right  circulation  about  the  qubit  loop,  and  its  net  magne¬ 
tization  is  discernible  using  a  dc  SQUID  [14].  As  a 
quantum  object,  the  potential  wells  exhibit  quantized 
energy  levels  corresponding  to  the  quantum  states  of 
the  macroscopic  circulating  current  [13,18].  These  levels 
may  be  coupled  using  microwave  radiation  [4,12],  and 
their  quantum  coherence  has  been  experimentally  dem¬ 
onstrated  [8]. 

Tuning  the  flux  bias  away  from  /  =  1/2  results  in  the 
asymmetric  double -well  potential  illustrated  in  Pig.  1(c). 
The  three  states  in  the  left  well  constitute  the  supercon¬ 
ductive  analog  to  the  atomic  A  system.  Using  tight- 
binding  models  with  experimental  PC  qubit  parameters 
[12-14]  at  a  flux  bias  /  =  0.5041,  we  estimate  the  inter¬ 
well  resonant-tunneling  rates  for  states  |1),  |2),  and  |3)  to 
berj"^^  ~  (1  ms)“*,r2'^'  ~  (1  /rs)“*,  and  ~  (1  ns)“\ 
respectively.  An  off-resonant  biasing  of  state  |3)  de¬ 
creases  its  interwell  tunneling  rate  to  order  (100  /us)“^ 
[23];  the  off-resonant  biasing  of  states  |1)  and  |2)  [as  in 
Pig.  1(c)]  will  also  significantly  decrease  their  interwell 
tunneling  rates.  In  addition,  the  intrawell  relaxation  rate 
at  a  similar  flux  bias  was  experimentally  determined  to 
be  =»  (25  /rs)“*  [12]  and,  presumably,  < 

P^mtra)  q’jjgj-efore,  the  “qubit  states”  |1)  and  |2)  are  effec¬ 
tively  metastable  with  respect  to  the  resonantly  biased 
“readout  state”  |3).  Since  Ej'^^  »  [12,13,18],  a 

particle  reaching  state  |3)  will  tunnel  quickly  to  state 
|4),  an  event  that  is  detectable  by  a  dc  SQUID. 
Alternatively,  for  slower  detection  schemes,  one  may 
detune  states  |3)  and  |4)  and  then  apply  a  resonant 
77  pulse  to  transfer  the  population  from  state  |3)  to  |4). 

Transitions  between  the  quantized  levels  are  driven  by 
resonant  microwave-frequency  magnetic  fields.  Assum¬ 
ing  fhe  Rabi  frequencies  be  much  smaller  fhan  all 

level  spacings  (U;|,  fhe  syslem-fleld  infer- 

acfion  may  be  wrillen  wifhin  fhe  rofafing  wave  approxi- 


mafion  (RWA)  [24], 


0 

0* 

iil2 

0* 

iil3 

12 

0 

f^23 

13 

^23 

-/E3 

(2) 


in  which  fhe  decay  from  sfafe  |3)  is  freafed  phenomeno¬ 
logically  as  a  non-Hermifian  mafrix  elemenf  [24,25].  Por 
small  microwave  perfurbafions  of  fhe  fruslralion.  A/,  fhe 
associafed  Rabi  frequencies  are  =  2TT\f{aEj/h)  X 
{p  \  sin(277/  -f  24>„)\q)\  numerical  simulafions  of  fhe  ma- 
frix  elemenfs  are  consisfenf  wifh  recenf  experimenfal 
resulfs  [see  capfion  Pig.  1(c)]  [12,13,18].  A  qubif  inifially 
in  |1)  can  be  prepared  in  a  superposifion  sfafe  = 
Ci|l)  +  C2|2)  by  femporarily  driving  fhe  012  held. 
Applying  fhe  013  (EI23)  then  allows  the  population 
of  state  1 1)  (|2))  to  be  read  out  through  a  transition  to  state 
1 3)  followed  by  a  rapid  escape  to  the  right  well  (a  readout 
scheme  also  utilized  by  single -junction  qubits  [7]). 

Alternatively,  one  may  achieve  S-EIT  in  a  supercon¬ 
ductive  A  system  that  is  prepared  in  state  Ifli'’)  by  simul¬ 
taneously  applying  the  microwave  fields  O13  and  O23 
such  that 


^13  ^ 

O23  Cl 


(3) 


Under  this  condition  (with  On  =  0),  the  state  is  an 
eigenstate  of  in  Eq.  (2)  with  eigenvalue  zero;  in 

this  “dark  state,”  the  SQC  becomes  transparent  to  the 
microwave  fields.  As  in  conventional  EIT,  the  amplitudes 
for  the  two  absorptive  transitions  into  |3)  have  equal  and 
opposite  probability  amplitudes,  leading  to  a  destructive 
quantum  interference  and  no  population  loss  through  the 
readout  state  |3).  Thus,  S-EIT  provides  a  means  to  con¬ 
firm,  without  disturbing  the  system,  that  one  had  indeed 
prepared  the  qubit  in  the  desired  state. 

In  a  practical  SQC,  there  will  be  imperfect  prepara¬ 
tion  as  well  as  decoherence  of  the  state  IflC),  and  this 
must  be  measured,  characterized,  and  minimized  for 
quantum  information  applications.  S-EIT  is  a  sensitive 
probe  for  this  purpose,  since  deviations  in  the  ampli¬ 
tude  and/or  relative  phase  of  the  complex  coefficients 
c,  from  the  condition  established  in  Eq.  (3)  result  in  a 

small  probability  |(cini3  +  C2^2'i)/^\^  (where  Ei  = 

^|fii3p  -f  IEI23P)  of  the  SQC  being  driven  into  the  read¬ 
out  state  |3)  on  a  time  scale  ~E3/n^.  In  general,  there  are 
two  categories  of  decoherence:  loss  and  dephasing.  Eoss 
refers  to  population  losses  from  the  metastable  states  |1) 
and  |2),  and  it  is  present  in  an  SQC  due  to  both  intrawell 
and  interwell  energy  relaxation.  Dephasing  refers  to  in¬ 
teractions  of  the  SQC  with  other  degrees  of  freedom  in 
the  system  that  cause  the  relative  phase  between  Cj  and  C2 
to  diffuse.  Both  types  of  decoherence  act  to  drive  even  a 
perfectly  prepared  state  |''E)  out  of  the  dark  state  defined 
by  Eq.  (3). 

We  describe  the  system  with  a  3  X  3  density  matrix, 
where  the  diagonal  elements  describe  the  populations. 
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and  Pij,  i  A  j  describe  the  coherences  between  levels.  In 
the  presence  of  the  EIT  fields  013  and  [I23  (with  012  = 
0),  the  Bloch  equations  govern  the  evolution  of  the  density 


matrix  [24]: 

i  i 

P  11  =  “TiPii  “2^13^31  +  2^i3Pi3>  (4) 

P22  =  “r2P22  “  2^23P32  +  2^23P23’  (5) 

P33  =  “r3P33  +  -O13P31  —  -O13P13  + 

i 

“2^23P23>  (6) 

P  12  =  “712P12  “  2  ^13^32  +  2^23P13>  C^) 

P  13  =  -ri3Pl3  +  ^f^t3(Pll  “  P33)  +  ^023Pi2,  (8) 

P23  =  “y23P23  +  2^23(P22  “  P33)  +  2^13P21-  (9) 


The  remaining  three  elements’  equations  are  determined 
by  p*j  =  pji-  The  decoherence  rates  =  (E,  +  ^ j)/'!  + 
.y(deph)  jQgg  dephasing  contributions.  We 

concentrate  on  the  regime  in  which  the  readout  state 
escape  rate  r3  =  (1  ns)“^  =  (27r)  130  MHz  dominates 
all  other  loss  and  dephasing  rates;  thus  713  ~  723  ~ 
r3/2.  Eurthermore,  we  assume  the  dephasing  rate 
yCdeph)  metastable  state  losses  Ej  and  E2,  set¬ 
ting  El  =  E2  =  0  and  712  ~  Theoretical  esti¬ 

mates  of  dephasing  rates,  such  as  7i2*^*’*’\  in  multilevel 
systems  were  recently  obtained  in  Ref  [26]. 

We  illustrate  an  S-EIT  decoherence  probe  example  by 
applying  EIT  fields  013  =  O23  =  (27r)  150  MHz  to  the 
dark  state  I^E)  =  (|1)  —  |2))/v^  and  numerically  inte¬ 
grating  Eqs.  (4)-(9).  With  712  =  0  the  system  is  station¬ 
ary  and  no  population  is  driven  into  |3)  (P33  =  0);  when 
we  include  a  dephasing  rate  712  =  (27r)  5  MHz,  P33  is 
small  but  nonzero  [Eig.  2(a)].  It  exhibits  a  rapid  initial 
rise  with  transitory  oscillations  [see  inset  Eig.  2(a)], 
reaching  its  maximum  value  p^™”^  within  about  = 
4  ns.  This  is  followed  by  a  smooth  decay  with  a  1/e  time 
of  about  80  ns.  The  solid  curve  in  Eig.  2(b)  traces  the  total 
population  P  =  p^  +  P22  +  P33  remaining  in  the  system 
as  a  function  of  time.  When  the  excited  state  maximum 
P3™’‘^  is  reached,  the  total  remaining  population  is 
JPiTss)  =  0.973.  In  contrast,  the  dashed  line  in  Eig.  2(b) 
illustrates  the  rapid  population  loss  expected  when  the 
same  fields  are  applied  to  the  state  1"^^)  =  (|1)  +  |2))/V2 
[tt  out  of  phase  with  the  dark  state  in  Eq.  (3)].  In  the 
absence  of  S-EIT  quantum  interference,  the  entire  popu¬ 
lation  is  lost  on  a  time  scale  ~  4  ns. 

We  now  use  Eqs.  (4) -(9)  to  show  how  measuring  the 
slow  population  loss  in  S-EIT  can  be  used  to  extract  the 
decoherence  rate  712.  The  elements  P33,  P13,  and  P23  in 


t(gs)  t(gs) 

FIG.  2.  (a)  Population  P33  versus  time  with  EIT  fields  H13  = 

fl23  =  (27r)  150  MHz,  an  initial  dark  state  pu  =  P22  =  0.5 
and  P12  =  —0.5,  and  with  a  dephasing  rate  712  = 
(2tt)  5  MHz.  The  inset  shows  P33  for  early  times.  P33  exhibits 
a  rapid  rise  to  a  plateau  after  a  short  time  followed  by  a 
much  slower  decay,  (b)  The  total  population  P  remaining  in  the 
system  versus  time  for  the  same  simulation  (solid  curve).  The 
dashed  curve  shows  the  population  for  the  out-of-phase  case 
Pii  =  P22  =  0.5  and  pi2  =  0.5  discussed  in  the  text. 


Eqs.  (6),  (8),  and  (9)  are  damped  at  a  rapid  rate  ~E3, 
allowing  their  adiabatic  elimination  [25,27];  we  solve  for 
their  quasisteady  state  values  by  setting  P33  =  P13  = 
P23  =  0.  This  approximation  is  accurate  once  initial  tran¬ 
sients  have  passed  and  the  plateau  value  P3™'  has  been 
reached.  Using  these  results  in  Eq.  (7)  yields  an  equation 
for  P12  with  a  strong  damping  term  YP'/Yj,,  and  it  too  can 
be  solved  for  its  quasisteady  state  value.  In  the  limit 
7l2E3/a2  1  we  obtain 

Pnit)  -  -  ^^^(1  -  ^^)[Pu(0  +  P22m.  (10) 

The  ratio  27i2E3/fi^  represents  the  small  fractional  de¬ 
viation  of  P12  from  its  dark  state  value.  There  is  a  com¬ 
petition  between  the  “preparation  rate”  YI^/Yt,  (which 
constantly  acts  to  drive  the  system  into  the  dark  state)  and 
the  decoherence  rate  712  (which  attempts  to  drive  it  back 
out). 

Plugging  our  adiabatic  solutions  for  pi3,  P23,  and 
Eq.  (10)  into  Eqs.  (4)  and  (5)  reveals  that  deviations 
from  the  dark  state  cause  a  loss  of  the  population  P  at  a 
rate  R  =  2712(0^31123/12"^)^.  Since  we  assumed  all 
population  is  lost  through  |3)  via  the  decay  term  P33E3 
in  Eq.  (6),  we  can  equate  these  two  rates.  When  the 
maximum  P3™^  is  reached  sufficiently  fast,  the  popula¬ 
tion  P  is  still  approximately  unity  and  this  yields 


(max) 

P33 


O13O23  7i2 

04  T7' 


(11) 


So  long  as  the  loss  during  the  initial  transient  time 
t  <  Ess  is  negligible,  the  population  will  follow  a  simple 
exponential  decay  P{t)  =  exp(  — p3™"‘’‘'E3t)  [as  in 
Eig.  2(b)],  and  the  dephasing  rate  712  can  be  easily 
extracted.  The  time  is  generally  the  greater  of  the 
preparation  time  ~E3/0^  and  the  inverse  of  the  decay 
rate  I/E3;  the  loss  up  to  f  =  will  be  ~p3““^E3Ess  ~ 
(20j3023/04)  max(7i2/E3,  7i2E3/0^),  requiring  both 
ratios  in  the  max(-  •  •)  argument  to  be  small.  The  first  ratio 
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FIG.  3.  (a)  The  maximum  plateau  value  ^3”“^  for  different 

yi2  (circles).  The  solid  curve  shows  the  prediction  (11).  (b)  The 
remaining  population  P  =  Pn  +  P22  +  P33  at  the  time  the 
plateau  is  reached  for  the  cases  in  (a). 


will  generally  be  small  for  parameters  of  interest  while 
the  second  can  be  made  small  by  choosing  an  appropriate 
field  strength  fl.  When  these  conditions  are  satisfied,  fhe 
mefhod  ouflined  here  can  be  used  fo  measure  712-  Nofe 
also  lhal  if  one  applies  fhe  S-EIT  fields  for  a  lime 
Ts-eit  ^  Yjs  and  observes  no  loss,  Ihen  fhe  prepared 
superposilion  I'P)  is  unperlurbed  and  one  has  slill  ob- 
lained  an  upper  bound  on  fhe  decoherence  rale  712- 

We  have  performed  a  series  of  numerical  simula- 
lions,  varying  712,  to  validate  the  preceding  approach. 
Eigure  3(a)  indicates  versus  712  and  compares  the 

numerical  results  with  the  analytic  estimate  [Eq.  (11)]. 
The  agreement  is  good  for  7j2  <  (27r)4  MHz,  which  cor¬ 
responds  to  27i2r3/fl^  <  0.056.  Higher  dephasing  rates 
compete  more  with  the  preparation  rate,  making  the 
adiabatic  elimination  approach  used  to  obtain  Eqs.  (10) 
and  (11)  less  valid.  In  such  cases,  one  observes  a  signifi¬ 
cant  loss  of  P  by  the  time  Pj™^^  is  reached,  as  illustrated 
in  Eig.  3(b). 

In  employing  the  RWA  in  Eqs.  (4) -(9),  we  have 
ignored  the  field-induced  couplings  of  far  off-resonanl 
Iransilions,  which  can  drive  fhe  system  oul  of  fhe  dark 
slate  [e.g.,  in  Eig.  1(c):  allhough  deluned by  |A|  ~  4  GHz, 
023  weakly  drives  1 1)  |2)].  We  have  performed  calcu- 

lalions  of  fhe  evolulion  including  all  such  couplings  for 
our  paramelers  and  found  Ihey  induce  an  addilional  ab- 
sorplion  info  |3)  al  a  rale  Roff-res  =  (271")  19  kHz.  Such 
Iransilions  will  nol  effecl  our  S-EIT  measuremenl  for 
qubil  decoherence  rates  712  ^  ^off-res-  The  effecl  is 
analogous  lo  ac  Slark  shifts  and  Rayleigh  scallering  in 
atoms  [28],  and  il  has  been  shown  [25]  lhal  Ihe  loss  rate 
scales  as  ~  r3(0^/A^),  0"^/r3A^,  which  inlro- 

duces  a  maximum  allowed  field  intensify  O^  for  a  given 
deluning  A  and  decoherence  rate  712- 

We  have  proposed  using  Ihe  superconduclive  analog  to 
EIT  (S-EIT)  to  demonslrale  macroscopic  quanlum  inter¬ 
ference  in  superconduclive  quanlum  circuils.  We  have 
shown  how  S-EIT  can  be  used  to  measure,  wilh  a  single 
pulse  of  Ihe  Iwo  S-EIT  fields,  whelher  a  parlicular  super¬ 
posilion  of  melaslable  energy  levels  (a  qubil)  has  been 
prepared.  The  technique  is  dislinguishable  from  previous 
slate  measuremenl  schemes  [7]  in  lhal  S-EIT  ideally  does 
nol  dislurb  Ihe  system,  preserving  ils  quanlum  coherence 


when  il  has  been  prepared  in  Ihe  desired  slate. 
Eurlhermore,  we  have  shown  how  S-EIT  can  very  sensi- 
lively  probe  small  qubil  errors  due  to  decoherence  or 
imperfecl  slate  preparalion,  and  we  have  oblained  ana- 
lylic  expressions  for  Ihe  field  slrenglhs  required  to  mea¬ 
sure  Ihe  qubil  dephasing  rate. 
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